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The X(3872) is universally accepted to be an exotic hadron. In this letter we assume that 
the A" (3872) is a D°D*° molecule, as claimed by many authors, and attempt an estimate ol its 
. prompt production cross section at Tevatron. A comparison with CDF data allows to draw rather 

^sL^ ■ compelling quantitative conclusions about this statement. 

^ \ PACS: 12.39.-x, 12.39.Mk, 13.75.-n 

^ n 

Introduction. Since the discovery of the X(3872) resonance by Belle [If it was soon realized that, despite its 

prominent decay mode were into J/i/jp, this particle could not have been identified as a standard charmonium 
00 1 excitation 0. The mass of the X being so close to the D°D*° threshold suggested to many authors that it 
could be a neat example of a hadron molecule 0] . The idea of molecules of hadrons dates back to the work by 
Fermi and Yang [J] in which the pion is interpreted as a nucleon-nucleon bound state. The case of the ^(3872), 
in this respect, is rather peculiar though: the D°D*° molecule constituting the X(3872) would be characterized 



■ by an extremely small binding energy. 

After Belle, CDF and DO confirmed the X(3872) in proton-antiproton collisions Later, BaBar Q found 

the X(3872) in the same B decay observed at Belle. It seems odd that a such a loosely bound molecule could 
be produced promptly (i.e. not from B decay) in a high energy hadron collision environment. 

This was one of the motivations to consider the possibility that the AT (3872) could be, instead of a molecule, 
a 'point-like' hadron resulting from the binding of a diquark and an antidiquark [g], following the interpretation 
proposed by Jaffe and Wilczek 0] of pentaquark baryons (antidiquark-antidiquark-quark). The drawback of 
the diquark picture is the proliferation of states predicted and the little insight on selection rules that could 
qq ' explain why many of these states are not yet seen (for a recent account on this aspect see e.g. |10|). 

In this letter we address the problem of estimating the prompt production cross section of the A' (3872) at 
the Tevatron and compare our results with the most recent observations by the CDF experiment allowing us 
to obtain a lower bound for it. We will assume that the AT(3872) is a D°D*° hadron molecule. For an earlier 

■ attempt in this direction [11 1 
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We employ standard Monte Carlo tools as Herwig [11] and Pythia [l3[ to compute the production of D and 
D* hadrons in proton-antiproton collisions at Tevatron energies. Open charm meson pairs produced are ordered 
as a function of their relative 3-momentum and of their center of mass momentum. We select those which pass 
the kinematical cuts used in the analysis made by the CDF collaboration. 

As we will explain below, we can estimate an upper bound for the theoretical cross section and a lower bound 
for the experimental one. The comparison of the two should give a qualitative answer whether the production 
of AT (3872) is exclusively due to the formation of a molecular bound-state. In the following considerations we 
will not make use of any particular model. The only model dependency in our calculations is that hardwired 
in the hadronization schemes of Pythia and Herwig. Indeed, to add weight to our conclusions, a comparison 
between the results obtained with both MC's will be carried out. 

CDF data on prompt A^(3872) production. The CDF collaboration soon after the observation of the 
A"(3872) in the J/i{)ir + ir~ channel measured the fraction of promptly produced AT(3872) to be 83.9±5.3% (l4| . 
In the same conference note we find the yields of i>(2S) — ► J/V>7r + 7r - and A^(3872) — > J/tpir + Tr~ candidates, 
using nearly identical selection criteria (indeed the selection differs by just an additional cut on the dipion 
invariant mass for the X(3872) case which is reported to produce little effect on the X(3872) signal), as well 
as the measurement of the prompt fraction of ip(2S) candidates. From these, and assuming the same detection 
efficiency for ip(2S) and AT(3872), which is presumably true within a factor of two but has never been reported 
as such by the CDF collaboration, we can roughly estimate, taking from [l5[ the B(-0(2S') — > J/ilnr + n~): 

<j(pp -> X(3872) + All) prompt x g(X(3872) -> J/V^+tt')) =i7±Qi „ (1) 
a(pp -> ^(25) + All) . . o ( ) 

where the uncertainty includes statistical and systematic uncertainties but do not attempt to gauge the assump- 
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tion of equal detection efficiency. The acceptance of the ip(2S) and A(3872) candidates is not specified in [Tj], 
but from the CDF II detector geometry and the indicated candidates selection we can conservatively assume 
that the above ratio applies for p± > 5 GeV and \y\ < 1. 

To derive an absolute A(3872) cross section, we can use the recently published jl6j ip(2S) Run II prompt 
differential cross-section measurement integrating from pt > 5 and taking B(-0(2S) — > ) from [15] : o~(pp — > 

V>(2S) + All) = 67 ±9 nb for p ± (ip(2S)) > 5 GeV, \y(^(2S))\ < 0.6. 

Assuming that both A(3872) and ip(2S) have the same rapidity distribution in the range \y\ < 1, we can 
finally estimate a lower bound on the prompt prodution cross section of A (3872) as: 

a(pp -> A(3872) + All) min > a(pp -> A(3872) + All) x 6(A(3872) -> J/^tt+O = 3.1 ± 0.7 nb (2) 

in p±(X) > 5 GeV, |y(A)| < 0.6. We aim at comparing c™™ with some upper bound on the theoretical 
determination of the cross section: er^ ax - 

Estimating an upper bound for <7 t h- Let us suppose that A (3872) is an S-wave bound state of two 
D mesons, namely a l/V2(D°D*° + D°D*°) molecule (we will use the shorthand notation D°D*°). Such a 
molecule has the correct I ++ quantum numbers of the A(3872). The A(3872) prompt production cross section 
at the Tevatron could be written as: 



a(pp -> A(3872)) 



d 6 k{X\DD* (k)) (DD* (k) \pp) 



2 



d 3 k(X\DD* (k)) (DD* (k) \pp) 

n 
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< [ d 3 k\^(k)\ 2 f d 3 k\(DD*{k)\pp)\ 2 < [ d 3 k\(DD*{k)\pp)\ 2 ~ a(pp-+ A(3872)) max (3) 
Jn Jn Jn 

where k is the relative 3-momentum between the D(pi), D*(y>i) mesons. tp(k) — (X\DD*(k)) is some nor- 
malized bound state wave function characterizing the A(3872). 1Z is the integration region where ^(k) is 
significantly different from zero. The first inequality is the Schwartz inequality: the e qual sign holds only when 
-0(k)/ (DD*ik)\pp) is equal to a constant almost everywhere in the integration region [Tj]. The matrix element 
(DD* (k)\pp) can be computed using standard matrix-element/hadronization Monte Carlo programs like Herwig 
and Pythia. To do so, we require our MC tools to generate 2^2 QCD events with some loose partonic cuts. 
Configurations with one gluon recoiling from a cc pair, are those configuration expected to produce two collinear 
charm quarks and in turn collinear open charm mesons. The parton shower algorithms in Herwig and Pythia 
treat properly these configurations at low p± whereas they are expected to be less important at higher p± . We 
will discuss also these processes at the parton level. 

In the hadron samples produced by the shower/hadronization algorithm we list the events containing D°D*° 
as a function of their center of mass relative momentum. At this level, the only cuts are those on partons: 

p part > g | y p art | < g j f 

more than one D°D*° pair is found in the event, we select the pair having the 

smaller relative 3-momentum k. 

We tune our MC tools on CDF data on D°D*~ pair production cross section distributions in the Acf> variable, 
(j> being the azimuthal angle in the transverse plane to the beam axis z; see Fig. [T] and [2j 

As for the determination of the region 1Z in ([3]) we estimate it having in mind a naive gaussian ansatz for 
the bound state wave function. It is straightforward to estimate the momentum spread of the gaussian by 
assuming a strong interaction Yukawa potential between the two D mesons. Given that the binding energy 
S is £ ~ M x - M D ~ Mjj* = -0.25 ± 0.40 MeV we find that r ~ 8 fm (8.6 ± 1.1 fm) and applying the 
(minimal) uncertainty principle relation, we get the gaussian momentum spread Ap ~ 12 MeV. One can check 
that changing the value of g 2 /47r ~ O(10) has a small effect on the size of 1Z. 

Given the very small binding energy we can estimate k to be as large as k ~ y / 2[i(~ 0.25 + 0.40) ~ 17 MeV 
or of the order of the center of mass momentum k = a/ X(m x , m 2 D , m*^)/2mx — 27 MeV. Given these consid- 
erations, we can restrict the integration region to a ball 1Z of radius 1 ~ [0, 35] MeV. 

Since we assume that D — D* interactions have a range of ~ l/m^, we expect a relative orbital angular 
momentum i < fc/m^, i.e., we can only allow S'-wave resonance scattering. Moreover we expect that D°D*° 
is a rather narrow object, its width being almost equal to the width of its D* component: T ~ 65 KeV. This 
is compatible with the determination of Belle and BaBar which find that the width of the A (3872) in the 
J/4>P channel is < 2.3 MeV at 90% C.L.. On the other hand attractive potentials do not generate such sharp 
resonances in S'-wave. In higher partial waves the centrifugal angular momentum barrier allows the formation of 



1 Which corresponds to a fco of the Gaussian at ~ 27 MeV and a spread of +12 MeV. 
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D" :: |y|<l :: 5.5<p i <20 GeV 
D'~ :: |y|<l :: 5.5<p x <20 GeV 
Herwig (p/" 1 > 2 GeV) 
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D° :: |y|<l :: 5.5<p^<20 GeV 
D*~ :: |y|<l :: 5.5<p i <20 GeV 
Pythia (pj& >2GeV) 



FIG. 1: The D D*~ pair cross section as func- 
tion of A(f> at CDF Run II. The transverse mo- 
mentum, p±, and rapidity, y, ranges are indi- 
cated. Data points with error bars, are compared 
to the leading order event generator Herwig. The 
cuts on parton generation are p^ rt > 2 GeV and 
\y p | < 6. We have checked that the dependency 
on these cuts is not significative. We find that we 
have to rescale the Herwig cross section values by 
a factor -KHerwig — 1.8 to best fit the data on open 
charm production. 



FIG. 2: The same as in Fig. [T]but using Pythia. 
We find that we have to rescale the Pythia cross 
sections by a factor Kp yt hin — 0.74 to best fit the 
data on open charm production. In both cases the 
agreement of the Monte Carlo distribution with 
data is remarkable. 



bound metastable states. Although as the D°D*° molecule has to be a 1 ++ state, we would need the first even 
parity wave, namely D-wave! Indeed in higher partial waves one can estimate the width of the resonance to be 
r ~ A£ ~ £o(ka) 2t ~ 1 , a being the range of the interaction. Then, in the D— wave case, there could have been a 
sharp resonance that we do not expect in the S— wave case where we necessarily are. Other molecule formation 
mechanisms under study, namely Feshbach resonances, could explain the narrowness of theses states (l8j . 

Results. As shown in Figs. [T] and [2 we can reproduce rather well the cross section distributions in azimuth 
intervals A0 for open charm production at CDF (see for example [l!| and the relative CDF internal notes), 
provided that we adopt some rescaling factors as to get the right normalizations. 

We have used Herwig and Pythia to compute hadron final states from 2-^2 QCD parton processes reaching a 
Monte Carlo luminosity C ~ 100 nb -1 . In Fig.[3]we show the integrated cross section as a function of the center 
of mass relative momentum in the D°D*° molecule obtained using Herwig. To get the minimal experimental 
value of a ~ 3.1 ± 0.7 nb we need to include D°D*° configurations having up to k IC \ = 205 ± 20 MeV. Molecule 
candidates in the ball of relative momenta 1Z can account only for 0.071 nb. Repeating the same calculation 
with Pythia, see Fig. [H we get fc re i = 130 ± 15 MeV whereas in 1Z we integrate 0.11 nb. 

Simulating the real experimental situation of prompt production of A (3872) at CDF would require a further 
increase of just a factor of 10 4 in the Monte Carlo luminosity which is extremely CPU demanding. Yet, in 
consideration of the stability of our results, we do not expect significant variations from what here observed. 

In conclusion we study gcc events with one gluon at p± > 5 GeV recoiling from the cc pair which in turn 
can hadronize into open charm mesons very close in phase space. We perform this computation at the parton 
level using ALPGEN [2(| and assuming that the fragmentation functions into open charm mesons to be set 
to one. This corresponds to an upper bound estimation. The results obtained point at a definitely negligible 
contribution from these configurations, being in the range of few picobarns. 

Conclusions. We have simulated the production of open charm mesons in high energy hadronic collisions 
at the Tevatron. The generated samples have been examined searching for D and D* mesons being in the 
conditions to form, through resonant scattering, bound states with binding energy as small as ~ 0.25 MeV. 
These A (3872) candidates have been required to pass the same kincmatical selection cuts used in the CDF 
data analysis. This allows to estimate an upper bound for the theoretical prompt production cross section of 
A (3872) at CDF. Averaging the results obtained with Pythia and Herwig we find this to be approximately 
0.085 nb in the most reasonable region of center of mass relative momenta [0,35] MeV of the open charm 
meson pair constituting the molecule. This value has to be compared with the lower bound on the experimental 
cross section, namely 3.1 ± 0.7 nb, extracted from CDF data The intuitive expectation that S— wave resonant 
scattering is unlikely to allow the formation of a loosely bound D°D*° molecule in high energy hadron collision 
is confirmed by this analysis. 



3 




0.0 



0.0 



0.2 



0.4 



0.6 



0.8 



LO 



0.0 



0.2 



0.4 



0.6 



O.S 



1.0 



* rel (GeV) 



FIG. 3: The integrated cross section obtained 
with Herwig as a function of the center of mass 
relative momentum of the mesons in the D°D*° 
molecule. This plot is obtained after the genera- 
tion of 55 x 10 9 events with parton cuts pP art > 
2 GeV and |y part | < 6. The cuts on the final D 
mesons are such that the molecule produced has 
a p± > 5 GeV and \y\ < 0.6. 



FIG. 4: Same plot as in Fig. [3] but using Pythia. 
We show these curves in a wide range of k IC \ to 
give an idea of the remarkable Monte Carlo sta- 
bility against fluctuations achieved on account of 
the very high statistics used. 
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